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SU fdARY 



U8in~ Holzer's aethod of frequency calculation, the natural 
frequencies for the first two modes of torsional vibration of the 
wing were determined for a representative conventional airplane 
(B24-C; in the custctaary manner, the fuselage being considered as 
a rigid body. Next, using a method developed by N. 0. Ifyklestad 
of the Guggenheim Aeronautics Laboratory at the California In- 
stitute of Technology, combined with Holzer's method, the natural 
frequencies for the same two modes of vibration wore again deter- 
mined, but with the fuselage this time being considered as flex- 
ible. 

A comparison of results of the two methods indicates that in 
considering the fuselage as being flexible, a decrease in the 
natural frequency of torsional vibration may be expected. For 
the particular airplane selected, this decrease amounted to 6.68^ 
for the first node of vibration and to 39.1% for the second. 

Bie investigation reported in this paper was entirely theo- 
retical and was performed during the 1945-1944 school year at the 
Guggenheim Aeronautics Laboratory at the California Institute of 
Technology, Basadena, California under the direction and super- 
vision of Dr. N. 0. Myklestad, research associate in aeronautics 
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at the Institute. 



INTRODUCTION 



In the design of aodern aircraft for hig,her and higher speeds, 
the designers are beooaing increasingly acre interested in the 
problems of flutter and vibration. One of these problens is that 
of the torsional vibration of the wings, which is dependent upon 
a nunber of factors, such as (Ij the aaas distribution both spcui- 
wise and cliordwise of the wings thenselves and of all units sun- 
ported either on then or within them, {Z) the torsional stiffness 
of the wings, (3^ the torsional noJiient annlied to the wings at the 
shifting center of oressure by the air loads, (4> the coupling be- 
tween the wings in bending and the wings in torsion, the tor- 
sional aoTient applied at the root of the wings by a flexing fuse- 
lage, nnd (6j the effect of conpressibility as looal velocities 
over the wing approach the velocity of sound. It is believed to 
be coTcaon oracticc in the aircmft industry generally to consider 
all but the last two of the factors enuner«.ted above. 

Ihis paper then Ims as its objective the qiamtitative deter- 
mination of the effect on a representative large airplane of the 
fifth factor emnerated above, namely, the effect of the torsional 
mcrnent apnlied at the root of the wings by a flexing fuselage on 
the natural frequencies of the wing torsional vibration. 
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DEFIHITIO- OF SYMBOLS 



- Fuselage Tiass concentrated at any station n, lbs. seconds 
squared per inch. 

rt - Nuaber of any wing or fuselage station. 

X - f.fass moment of inertia of wing in. inch lbs. seconds squared 
or bending moment of inertia of fuselage in inches to the 
fourth power. 

- Convenient reference value of bonding moment of inertia for 
fuselage as a whole. 

E - 'Modulus of elasticity of fuselage bending material in lbs. 
per square inch. 

“ Angular deflection of wing at any station n in radians. 

tA - Frequency of vibration in radians per second. 

- 5^el length of wing or fuselage between stations n and 
n •+• 1 in inches. 

5„ - Shear at fuselage station n in lbs. 

- Bending moment at fuselage station n in inch lbs. 

Mb - ” " of fuselage tail at elastic axis. 

M; - ” t, r .. I. t. 1. ^ 

oCfl - slope of fuselage axis at any station n. 

" " ” tail at elastic axis. 

oCt - " " " nose " " " . 

- Deflection of fuselage at any station n in inches. 

" " " tail et elastic axis. 
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- Deflection of fuselage nose at elastic axis. 

- Change in slope from n to vs + i due to a unit force at vs. 

Vm* " "" " n"vn-*l"""" monent at n 

dp - " " deflection froa n to n + l due to a unit force at n. 

dM„ “ " " " + jnoaent at n 

- Slope of fuselage axis at extresio end of tsiil or nose. 

- Coefficients appearing in equation for O-r, ■ 



I 






If 



n » 



n 




- Total coupling moment introduced into wing by fuselage at the 
elastic axis. 

= I^j/To ” '^Jon-di-nensional symbol for fanelage bonding moment of 
inertia at any station n. 

bn=a„^,-a^= rtjiiJ— — iiii- - Non-dimensional symbol for increase in 

fuselage bending moment of inertia froi 
station n to station n + i . 

J - Variable distance from station n to any point in panel 
(between n and n+l J. 






K„ = X LI.- X K.l 

A - \ S : I 

, A-r\-i ^ 

g: = X G,i 

i=« 



U - Torsional 

n 



rigidity of wing in lb. 



inches per radian. 
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A»'ALY';iS 



3eoausa of the nature of this investigation and the diffioultiea 
involved in the aeasurencnt of the effect of one factor at a time 
on the torsional vibration of a wing, no experimental work was under- 
tedcen. Instead, the authors approached the problem frcm a ourely 
theoretical viewpoint, end the investigation was performed entirely 
on that basis. 

After a representative airplane for the investigation had been 
selected, it was necessary first to obtain the following information 
oonceniing it : 

A. The wing (Table I;, considering its mass anl the masses of 
all bodies either attached to it or stored within it as 
being concentrated at a mraber of stations along its span: 

(1^ the distance of each station from the wing root in 
inches, 

(2j the mass polar aionent of inertia about the elastic 
axis of the wing of the mass considered to be concen- 
trated at each station in lb-inches seconds squared, 
and 

(3 ; the rigidity in lb-inches per radian, or its 

reciprocal, of the wing in torsion between each eta- 
tion. 

8. The fuselage (Table II;, considering its mass and the masses 
of all- bodies either attached to it or stored within it as 
being concentrated at a n’jnber of stations along its length: 
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( 1 ^ the distance of oach station frtxa the fuselage nose 
in Inches, 

( 2 > the bending moment of inertia about a horizontal axis 
perpendicular to the longitudinal axis of the fuselage 
In in inches to the fourth power, end 
(3j the total ness considered as concentrated at each 
station in lbs. seconds squared per inch. 

After receiot of the required information for the wing, it wae possible 
to calculate the natural frequencies of the wing in torsion for as 
many modes of vibration as were desired, considering the wings as being 
built in to a stiff fuselage with an extremely high aouent of inertia 
compared with that of each station along the wing. This calculation 
was actually carried out for two modes of vibration following Holzer's 
■nethod as outlined cm pages 228 and 229 of Wef. 1, on example of which 
has been aopended to this paper as Table III with an explanation in- 
cluded in the aopendix. The results of this calculation have been 
tabulated in Table 'IV and plotted on Fig. 1, and show that the natuiral 
frequencies for the first two modes as detemined by this oaloulation 
are 33.67 and 71.70 radians per second res'ieotively. 

Ihis completed the first ohase of the investigatior 5 and with the 
required information for the fuselage then at hand, it was possible 
to proceed with the second phase, namely, the calculation of the natu- 
ral frequencies of the wing in torsion for as .many modes of vibration 
as were desired, considering the torsional moment applied at the root 
of tile wings by a flexing fuselage. Hie nroblesas immediately con- 
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frontini' the nuthora in this phuao of the investi£.ation were those 
of deteroining (ly the torsional aotaent oroduced at the root of a 
wing by a flexing fuselage and (2} the method of coupling this aoraent 
into the wing at its root. 

For the solution of the first of these probl«as a method devel- 
oped in Ref. 2 for the antisjacaetric bending of wings was applied 
to the flexing fuselage, considering the fuselage to be made up of 
two independent beams extending in opposite directions from the lo- 
cation of the clastic axis at the root of the wing. This aethod has 
the advantage of yielding Lmaediately the bending moment at apy par- 
ticular station along a cantilever beam and the slope of the beam at 
that station as linear functions of the normal displacement of the 
beam. Consequently, the procedure followed wae, first, to calculate 
the bending aocments and the slopes, at the location of the elastic 
axis at the root of the wing, of both the portion of the fuselage 
aft of this location and the portion of the fuselage forward of this 
location resulting from a unit downward displacement of the extreme 
end of both the tail and the nose. The bending moment at the elastic 
axis end the slope at that location of the after portion of the fuse- 
lage were designated as Mj, and ocj, respectively, euid of the forward 
portion of the fuselage as and oc't, respectively. 

Hext, since the fuselage is actually a continuous structure 
throughout its length, its slope on either side of the elastic axis 
must equal the slope on the other side of the elastic axis. This 
leads to the result that, since the initial diaplacenonts of both the 
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tail and the nose were taken tc oe. so-aibive dowrr/rards, in order for 
the slope forward of the elastic axis to equal that aft of the elastic 
axis, the slcpo forward of the elastic axis °<b laust be .nultiplied by 
the ratio This sane result would have been obtained had 

the initial diaplacenent of the nose been multiplied by this ratio 

• and since the bending moment developed is a linear func- 
tion of the displacement of the free end, the ben ling moment produced 

/ 

at the elastic axis b, the forward portion of the fuselage should 
also be multiplied by this same ratio. ”'e then have that, for the 
continuous fuselage, the bending moments at the elastic axis due to 
the after and forward portions of the ^uselago are given by the ex- 
nressions, Mj, and (-oc./oc,). Mb respectively. HoiYeTer, these two 
cenponents oppose one another; consequently, in order to detemine the 
total bending moiont frera the fuselage to be coupled into the root 
of the wing, one must be subtracted fren tho other. If the direction 
of Mt, is taken to be the positive direction, it nay readily bo seen 
then that = (- ocb/oc;)Mb= And since the bend- 

ing moment fron the fiiselage at the elastic axis enters the wing as 
a torsional monent, torsional moment produced at the root 

of the wing, the amount of it entering each side of the wing being 
assutting symmetrical twisting of the wing* 

?or the method of coupling this aouent into the wing at its root, 
one side of the wing was considered as a free body in tor'^ion with this 
torsional laaaert of i:M^ aoplied at its root. Assuming arbitrary unit 
angular deflections of the wlnj tip, the Holser’s calculations nado 
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during the first phase of this iavestigution yielded the following in- 
fornation for each frequency selected: 



(ly The bojrsional laoncnt derelopod vrf.thin the wing at the first 
station outboard of the root due to the rotational inertia 



forces within the wing ( froa Table III;, and 

(2; Tne angle of twist developed at the root of the wing ( 
froa Table III;. 

>inco in this calculation the torsional mofaent developed at any 
particular station and the angle- of tv/ist at that station are given 
as linear functions of the arbitrary angular deflection of the wing 
tip, any desired angle of twist at the wing root can be obtained by 
properly adjusting the arbitrary angular deflection of the wing tin. 
;ince the wing can be considered to be built-in to the fuselage, its 
angle of twist at the root should equal the slope of the fuselage at 
the wing’s elastic axis, and in order to obtain this an/le of twist 
at the root it is necessary to nultiply the original arbitrary angular 
deflection of the wing tip by the 'ratio davi.ng aaltiplied 

the original arbitrary angular deflection of the wing tip by this 
ratio, it is then necessary to multiply t’ne torsional moment devel- 
oped 7 n.thin tVie wing at the first station outboard of the wing root 
by this ratio also. Hence, this moment is then found to equal 



, aad adding tlus to the torsional aonont applied at 
the root of the wing by the fuselage, a re-idual torsional moment or 



torque is hhen the additional applied mcnent required in order tc force 





is found. This residual 
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t 



the wing to vibrate at the assuned frequency. Tnis value of the re- 
sidual torque is then plotted against the asauTied frequency; an’ the 
process is then repeated for other assuaed frequencies, one point 
on the plot being obtained for each assuned frequency. A ooaplete 
exanple of the calculation by this Tietho<^ for one assuned frequency 
has been appended to this paper as Tableo V(aj, V(b^, Vl(ay, VlCby 
and VII, with a brief explanation of then included in the appendix. 

The results of this calculation have been tabulated in Table VIII. 

After a sufficient number of points have been obtained, a curve 
may be drawn through then as has been done in Fig. 1. Again the 
points at which this curve crosses the frequency axis determine the 
natural frequencies of torsional vibration for the wing, for at these 
points the residual torque becemes zero, and hence the additional 
applied moment required to force the wing to vibrate at that frequen- 
cy also becomes zero. From Fig. 1 it can be seen t-^at the natural 
frequencies for the first two modes as determined by this calculation 
are 51.42 and 45.65 radians per secemd respectively, (lee Table IX 
for tabulation of final results. >• These aro reductions of 6.68^ and 
59.1%) respectively fron the frequencies of the first two modes found 
in the first phase of this investigation; Accordingly, it may be con- 
cluded that, whereas the consideration of a flexing fuselage has a 
small but appreciable effect on the frequency of the first mode of 
torsional vibration of the wing, it has a very decided effect in 
lowering the frequency of the second mode. Probably, considering the 
trend of the curves of Fig. 1, this same effect is carried on in pro- 
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^ression to subsequent -nodes of vibration; hence it is the studied 
opinion of the authors that this effect should be considered in the 
calculation of the natural torsional frequencies of the wing. 

Hie closing jxhase of this investigation was the determination of 
the fxiselage deflection curves for each of the -two modes of vibration 
deteimined above. This was accomplished with facility from the cal- 
culations involved in the determination of the bending aoaents and 
slopes of the forward and the after portions of the fuselage in the 
second phas* of this investigation. The deflection at any station of 
the fuselage is designated as , and oolimns so headed may be 
found in both Tables ''r (a; and \’l(by. Again, the values of 
given in Table VI (ay must be .aultiplied Oj the ratio in 

order to give them the correct magnitude with respect to those given 
in Tabic VI (b/. fuselage deflections are tabulated in Table X. A 
plot of the values of calculated for ‘he two natural frequencies 
found in the second phase of this investigation was made and has bren 
appended to this report as Fig. 2. A perusal of this figure will in- 
dicate thet tire deflection curves for the fuselage for the -two nodes 

t 

of wing torsional vibration aiv very similar, there being no reflex 
curvatures along the fuselage length in either case. Tn the first 
mode the nose deflection is about .one-seventh that of the tail whereas 
in the second mode it is almost twice that of the tail, fran which it 
can be seen quite readily that for a ri\en deflection of the nose the 
fuselage curvature will be much greater for the first mode than it 
will for the second. 
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Tno relative angular deflections of the wing at each station 
along its span can be determined very readily by referring to the col- 
umns headed /3 in Table III for the calculation for a rigid fuselage 
and in Table VII for the calculation for a flexing fuselage. These 
values must be multiplied by the ratio for each frequency 

selected in the calculation for a flexing fuselage, as has already 
been done in the detenoination of the residual torque acting on the 
wing, in order to detenaine the actual magnitudes corresponding to 
a unit dovnrvard deflection of the tail. This has been done and the 
results for the two modes tabulated in Table XI and plotted in Fig. 8. 

Fig. 6 is a schematic illustration of the two modes of vibra- 
tion, assuaing a unit downward deflection of the extreme tail in each 
case. 
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CCNCLUSI0S3 



In the case of the airplane investigated herein, the consideration 
of a flexing fuseleige has a saall but appreciable effect on the fre- 
quency of the first mode of torsional vibration of the wing, but it 
has a vejy decided effect in lowering the frequency of the second mode. 

The deflection curves for the fuselage for the first two modes 
of wing torsional vibration are very similar, there being no reflex 
curvatures along the fuselage length in either case. However, for a 
given deflection of the nose, the, fuselage curvature will be much 
greater for the first mode than it will be for the second. 

It must be understood that the above ocsiclusions apply only to 
the particular airplane which has been investigated herein. Fnis 
paper is not submitted with the intent to show that effects of similar 
magnitude can be expected for all airplanes, but simply that the effect 
should be investigated with the thought in mind that it might orove 
appreciable, oartioularly in the case of higher nodes. 



% 
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A^PLIIDIX 



I. calcui'vTion for the rigid FTRF.LAGE 

The method u^ed is outlined on pages 228 and 229 of Ref. 1, and 
is known as ’^folser'o method. The wing data (.'ee Table I.; furnished 
for the airplane in question assiaaed the mass moments of inertia X„ 
of the wing to be concentrated at 7 s-^anwise stations, the first and 
last stations being located at the tip and the root respectively. 

The notation used is the same as that for the fuselage and is demon- 
strated in Fig. 4. 

A positive (climbing; pitching angle at the tip ( n-i j of 
one radian was assumed, ( ^ = ; and for a given frequency, the 

inertia torque was calculated for station n = i . This 

inertia toroue multiplied by the torsional flexibility for panel 
length i, gave the amount of twist or 4he reduction in angle be- 
tween stations n=t and n = 2 . This angle of twist was then sub- 

tracted from to give , the angular deflection at station 
n = 2 . Knowing , the inertia torque at station n=2 , 

was then calculated. 

Ihe remainder of the table was completed in like manner until 

0-7 ^ 

the residual inertia torque wing root,) was found. 

Ihis value of inertia torque was tabulated in Table IV and plotted 
against cO in Pig. 1. The residual inertia torque is the shaking 
moment which would be required at the wing root to cause the wing to 
vibrate torsiooally at the assumed frequency a) . •..f c:. -e 
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of a rigid fuseleige and assuniing syiuiotrio torsional vibration, 
aust equal zero at e natural torsional frequency of the wing. Consequent- 
ly these natural frequencies can b« found by plotting ^ against 

cd to deternine points of intersection with the tO axis as was done 
in F’ig* 1. This calculation was carried out for the first two nodes, 
the final results aopearing in Table IX under "Higid Fuselage”. 

A sanple calculation for appears a* Table III. 



II. CALCULATION FOH THE FLEXIBLE FU5ELACfE 
(a,^ Outline of Method Usedi 

The method used here for the fuselage is that derived in Pef. 2 
for the antisymmetric bending of airplane wings. A brief resume of 
the method follows herewith. 

Using the notation demonstrated in Fig. 4, the shear at any sta- 
tion n is given by 

s. - M 

and the bending aoment at any station n is given by 

( 2 ; 

The slope at any station n+ l is given ly 
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arid the deflection at any station n + \ is given by 




(4y 



where 

= change in slope froa n+v to lo due to a unit force at n. 

change in slope from n + i to n due to a unit nanent at n. 

= char.ge in deflection from n+i to n iue to a unit force at n. 

c1m„ = change in deflection froa n+i to n duo to a unit nouent at n. 

Tire ncthod of calculation of these naraacters is outlined in sec- 
tion II (d/ of thi* appendix. 

.Substituting the expreesiesas for S„ and M^, frota equations (1; 
and (C; into equations (Sy and (a; 




/ 




At the end of the fuselage, ass me 



ot, - 



so that c>t^ = 



(Ray 




Continuing \rith eqiaatioas (.Ry and (C j in like manner yields 






(7; 
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isgt,. 






m 



« 



t 



I 














( 8 ; 



wh«r« coefficients ®’*’® independent of , 

one caapletc set beinc obtained for each frequency. Tlie method of 
determination of these coefficients is outlined in section II (b^ of 
this ao-'endix. For the present, these corfficients are assmed to 
bo known for any particular frequency tO 

For antisyaaetric bending, the case where the fuselage is being 
shaken by a shaking mocaent M cos a)t about the elastic axis of the 
wing, the deflection at the elastic axis is zero O and from 
equation (3/ <^ = • 7«ith this val le of 4* » deflec- 

tions may be found by means of equation (G;, as can the bending ‘ 
Dicment at any station n , 2 . 

In the Particular calculation ?dth which t'.iis paper is concerned, 
the two quantities desired are the bending moment caaing in fror.a the 
tail or nose , and the slooe of the fuselagp at the elastic axis 









o; 


= m. (x, -xj - 2 m, k, (x, - X J 


(lo; 


Putting 


(11; 


and o)*" = K;, 


(12; 


also ^ Z Xj 

5 - A. 


(13; 
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men ^ K] - ^ 2 43 (U; 

XS( 5** t ^ 

But 4] = Q(4-t t-t-^-,) 4 s.(4 + 4+-'4-.) +--&b.. ife., 



“■^,4 + 0 + 4^*^i ■^*^1'^ -^t>-iC^' 

A«bH _ 4. _ 

= .2U^>J 



" ^fc-3 



(1^; 



and sinilarly ~ 'P ^ ^-4 ^ 






. - 2[X.|&J -<^IU|K,1 



, = l ' S-i s- 

^'* - i»-A 



50 



M 



(16; 



T?cf erring to Table VI (bj- 



^ Ks is given by colonn (2/ and 



is given by colunn (6;, the second smEiatione ocoui^ 

s* » 

ring In coluris (3) and (7^ which colurans give 



I i«t-l S*r4 

■K,= ^[X.ZKJ 



and 



/ >»sV»-l S»A 

2[1. ZG-.] 

I 5*1 



(17; 

(18; 



Fron this it is seen that 

K - ^'b - Kb ^ (19; 

and froTi equation {7 J o^b= - 4 (2o> 

where and 4 given on line Q under coluans (^) and (8; 

respectively. 

The .method was repeated for the nose using the sane frequency 
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(Table to obtaixi and o-l , and the four ■'aluea thus d«- 

teraincd arere tabulated at the top of Table VII. 

..'ith an assuned poaitive (doivnn'ard^ deflection of one inch at the 
tail the total noraent introduced at the elastic onia b, the fuselage 
is given by 

Me = (o^b/ocL)Mb + Mb 

and the slope of the fuselage at the elastic axis is <oc^ . 



(by Oeterniuation of Coeff icients 

hith the original assumptions at the end of tho fuselage ^ 

and i , frail equations (7j and (0> 

ot^-= - fr. 



it is obriou. that 

I -f, = o - 1 k. = o , and 

fran equations (5a^ and (6a^ 

o<^ ^ ^ 

are obtained -6^= i and • 

dubstituting equations (7 j and (S; into equations {^:j and (6;, 






( 22 ; 



3^ oq'1-.tlnj^, tar.:' co;\tainii ^ fS and t’^'ose not ccataining on the 

fcy.'o sides of cnc’i of eqintions (21; and (22;, the following eqaat ons 
art obtained: 

ITq'iatlhg coefficieata of "p : 



L„ “ £ nn.k, Cx*-X„) (2S; 

l^n., = 2 m.k - m. k (x^- x„) (2-»; 



Equating consta-nt temn : 



L. - i 



+ a) Vp 















in., = fn + KL, ' Z 



(25; 
(20 ; 



Using substitutions and relations developed in (11;, (12;, (13;, (15; 
and ( It ; ; 



A -O-* S = A. 



+ Vf.zt^. + V 'Z, [i 5 ^.] 

-k. .t.L, -d.„ZK. -d„;i:[i.‘f K,-] 

A-' X - ' I 

A-W i-h-< ^ S»1 

t = f„ ■+ 2 4 2 u Z 

in., ^ in + + . 1 



(27; 

(28; 

(29; 

(30; 



/ill the coefficients , k,, » > f^nd can be found by pro- 

gressive calculation with the aid of a table (such as Tables VI (a; and 
VI (b;; based on those equations. 



(c; Procedure Used in filling out Tables Vl(a; and Vl(b; 



Paraneters Vf^ , , and calculated (see section 

Il(d; of this appendix; and written in the spaces indicated. Like- 
wise the values for were entered in the tables. 

Then, for a given frequent, values of were calculated and 

entered in the appropriate colunn. 

’?e:tt, in line n=i , the following values were entered in 
colufiins (i;, ( 4 ; and (S; respectively: 

k, * o ■^ 1 = ' ‘3'' “ ' 

and in line n = z under colinm (1^: kx= -?■, 

The table was then worked across fren left to right starting with 
line n-i then proceeding with line n *2 etc., each step being in- 
dicated in the colurm heading. 

In ermputing values to enter in col’-wms (1,., (3;, end ( 7 ), caie 
Taust remember to u.se inl'omation appearing in the preceed.ing line. 

The remainder of the steps are self explanatory, the desired 
quantities of the calculation being and o;^ . 

(d; Calculation of Para-ueterr. , and 

The fuselage data received (Table II. ^ indicated bending moments 
of inertia equal to eero at each end of the fuselage, but in order 
to more nearly approximate the probable moment of inertia distribution 
In the regions from the extreme ends to the next stations inboard, 
a tranesoidal distribution over these regions was assumed in botn 
oases with the end ox^inates approximately half the value of the next 
ordinates inboard. The assraed values were 

•I. (TAIL) = 500 in.^ I. (NOSE) = 3,500 in'^ 



- 2S 



A fuselage stiffness curve (such as fig. could b© dravm for 
the given and assuaed (endy values of bfading aoaent of inertia, where 
the bending moment of inertia is plotted against fuselage distance 
as abscissa, such that the areas between succeeding stations would be 
trapezoids. If is talren as a convenient refe'^vmce value of the 
bending moment of inertia for the fuselage as a whole, then at any 
point between stations m and n-t-i 



the moment area method: 

0 

^ ) - a.bj- 



(Sly 

(Z2) 

( 33 ; 



Referring to Tables V(ay and (by, EX^ was taken as io'° • Colixons 
(Ij and {e:J were filled in with values of end a„ from the data 
furnished. The value (a^fku,) in colusnn (5^ is determined from 



a„+, = + b„ 

Due to lack of availability of a seven nlace table of natural 
logarithms,- comraon logarithms were used and values converted to 
natural logarithms in coluons (6J and (9J by the relation 

•H- 

The remainder of the table is self explanatory and follows from 



- 26 - 



equations (51/, (52 J and (55 The desired nara.aetars appear in col- 
(0^, {12 J and (18j. 



{e J ■'rocedure ”3«d in Filling out Tabic VII. 

/gain using the sane frequency as was used in Tables VI (a j and 
Vl(by, Holzor’s calculation (as explained in section I of this apocn- 
dix; was repeated in Table VII, the desired quantities being the 
angular deflection of the wing at the root /3^ and the inertia torque 



at the next station outboard fjxra the root . This inertia 

torq\;e ndjustO ' sc as to the on^rular deflection oi th< »d.ng at 

the root equal to the slope of the fuselare at the elastic axis is 



the >y=»nd5n^;^ mrrient introduced at the elastic oris \ry the fuselage 
i: gives the residual torque actin' at the w:ng root whic!* would 

require^ to TiaJro the wing vibrato torsicnally at the chosen frequency. 



At a natural frequency of the systea tnis residual torque is zero. 
Table VIII gives a tabulation of th^ results of this calculation, 
i ig. 1 shows a nlot of this residual torque against cO and the natural 
frequencies ffirst and second nodes y occur when this curve crosses the 
oO axis. The final results are tabulated in* Table IX under '’Flexible 
} uselage". 




, and this adjusted inr-tia torque added to half 




?7 - 



HMT’cr.NCIS 



Befercnc© 1 
Reference 2 



Den ^arto^; "'tcohanical Vibrations”, McGraw-Hill 
Book Company, Inc.| Sew York and London, 1940. 

Jf. C. Hykelstadi "A Now Hethod of Calculatinc Kat- 
nre.l Modes of Uncoupled Bendin{^ Vibration of Air- 
plane Things and Other Types of Beans”. Journal 
of the Aeronautical "'•ciences, Vol. 11, Ko. 2 , 
April, 1D44. 
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~ \AQ.l 
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+ 0. 4 1 3 
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3-4e 
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1 1 C.S 


0.99178 1 
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3 
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1 .14-1 
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4 
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0.9GG309 


1.128 
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+ 0. 38 7 


s 


a 17 
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0.^42 
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- 0.00 1 


G 


1 33 
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-0 247594 


- 0. 1 02. 


1 


0 
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- 0. 03 1 


•• 0.49G930 


- 0.2 05 
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FIG. 1 - RESIDUAL TORQUE VS. FREQUENCY FOR RIGID AND fLEXIBLE FUSELAGE 
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FIG. 2 - F'JIELAG:-. DEFLErXIO'N CURVES 
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flG. 4 - SKETCHES SHaTING NOTATION USED IN BENDI'.G CALCULATIONS 




PISTANCr FROM elastic AXU 



FIG. 5 - SAMPLE FUSELAGE STIFFNESS CUr^F. A17D NOTATION USED 
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